Specific viscosity of neutron-rich nuclear matter from the relaxation time approach 
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The specific viscosity of neutron-rich nuclear matter is studied from the relaxation time approach 
using an isospin- and momentum-dependent interaction and the nucleon-nucleon cross sections taken 
as those from the experimental data modified by the in-medium effective masses as used in the 
IBUU transport model calculations. The relaxation time of neutrons is larger while that of protons 
is smaller in neutron-rich nuclear matter compared with that in symmetric nuclear matter, and this 
leads to a larger specific viscosity in neutron-rich nuclear matter. In addition, the specific viscosity 
decreases with increasing temperature because of more frequent collisions and weaker Pauli blocking 
effect at higher temperatures. At lower temperatures the specific viscosity increases with increasing 
density due to the Pauli blocking effect, while at higher temperatures it slightly decreases with 
increasing density as a result of smaller in-medium effective masses at higher densities. 

PACS numbers: 21.65.-f, 21.30.Fe, 51.20.+d 



I. INTRODUCTION 

The properties of nuclear matter under extreme con- 
ditions is one of the major problems in nuclear physics. 
Transport models are now widely used in studying the 
equation of states (EOS) of the nuclear matter formed in 
intermediate-energy heavy- ion collisions [l[ ■ It was found 
that the transverse flow from the experimental data can 
be reproduced by the transport model calculation us- 
ing a stiff EOS together with a momentum-independent 
potential or a soft EOS together with a momentum- 
dependent potential 0]. However, the extracted incom- 
pressibility from Giant Monopolc Resonance experiments 
is 231 ± 5 MeV indicating a soft EOS and the impor- 
tance of the momentum-dependent mean field potential. 
Recently, an isospin- and momentum-dependent inter- 
action together with the Isospin-dependent Boltzmann- 
Uehling-Uhlenbeck (IBUU) transport model has been ex- 
tensively used to study the dynamics in intermediate- 
energy heavy-ion collisions. Especially, extensive stud- 
ies have been done to extract the information on the 
isospin-dependent part of the EOS, i.e., symmetry en- 
ergy, from isospin diffusion Q, neutron/proton [f| and 
triton/ 3 He @ ratios and differential flows, and 7r~ /tt + ra- 
tio • For a recent review, I refer the readers to Ref. Q . 

Although hydrodynamic models have seldom been 
used in intermediate-energy heavy- ion collisions @ , they 
are widely used in relativistic heavy-ion collisions, where 
a very strong interacting matter called quark-gluon 
plasma (QGP) is formed. It has been found that this 
matter has a very small viscosity and behaves like a 
nearly perfect fluid. Below the temperature T c = 170 ~ 
180 MeV hadronization happens and the viscosity of the 
system increases. Although there have already been ex- 
tensive studies on the shear viscosity of the QGP fl0l - fl3T | 
and that of the relativistic hadron gas flil- 17[. only a 
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few studies have been devoted to the viscosity of nuclear 
matter formed in intermediate-energy heavy-ion colli- 
sions (T8T - [22l |. Even fewer studies have been devoted to 
the isospin dependence of the viscosity [23| . It is thus of 
great interesting to study the viscosity of the neutron-rich 
nuclear matter, especially its isospin dependence, formed 
in the intermediate-energy heavy-ion collisions as in the 
IBUU transport model. 

Many methods in the literature have been used to 
study the viscosity [12, [IH, 0, , among which directly 
using Green-Kubo formulas is the first-principle way of 
the study [2|| . In the present study the relaxation time 
approach is used [26l \2l\ . This approach is helpful in 
studying how the Pauli blocking, in-medium cross sec- 
tions and in-medium effective masses will affect the shear 
viscosity and can give qualitative ideas how the shear 
viscosity changes with density, temperature and isospin 
asymmetry of the nuclear matter. 

This paper is organized as follows. The isospin- and 
momentum-dependent interaction (here after 'MDI') is 
briefly reviewed in Sec. |TTJ The method to calculate the 
shear viscosity from the relaxation time approach is given 
in Sec. Mil In Sec. IIVI the results of the relaxation time 
and the specific viscosity are shown and the effects from 
isospin, temperature and Pauli blocking are discussed. A 
summary is given in SecfVl 



II. THE ISOSPIN- AND 
MOMENTUM-DEPENDENT INTERACTION 



The MDI interaction is a modified finite-range Gogny- 
like interaction [28j], and it was recently found that 
the effective nucleon-nucleon potential is composed of a 
zero-r ang e many-body term and a finite-range Yukawa 
term [23 ]. In the mean- field approximation, the single- 
particle potential for a nuclcon with momentum p and 
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isospin t in a nuclear matter medium is written as 
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In the above r = 1/2 ( — 1/2) for neutrons (protons); 
S = (p n — Pp)/p is the isospin asymmetry with p n (Pp) 
being the neutron (proton) density and p = p n + p p being 
the total density of the medium; f T {p) is the local phase 
space distribution function and it can be written as 



U(p) = dn T (p) 



exp 



(& t +Ur(p)-Pr)/T 
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where d = 2 is the spin degeneracy and n T (p) is the 
occupation probability, m is the nucleon mass, T is the 
temperature and p T is the proton or neutron chemical 
potential and can be determined from 



Pt = J fr(p) 



d 3 p 
(27p 
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From a self-consistent iteration method [30j , both the 
single-particle potential and the distribution function can 
be calculated numerically. The entropy density of the 
nuclear matter can thus be calculated from 



= — d / [n T lnn r + (1 — n T ) ln(l — n T )\ 



d 3 p 
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The detailed values of the parameters A u (x), Ai(x), 
a, B, C TiT , C T - r and A can be found in Ref. [28j and 
they are assumed to be independent of the temperature. 
These parameters lead to the binding energy E = —16 
MeV, the incompressibility Kq = 211 and the symme- 
try energy E sym « 31 MeV at the saturation density 
pa = 0.16 fm~ 3 . The parameter x is used to mimic 
the density dependence of the symmetry energy while 
keeping the properties of symmetric nuclear matter un- 
changed. Comparison between the results from the IBUU 
transport model calculation with the experimental re- 
sults of isospin diffusion data leads to — 1 < x < at 
subsaturation densities 0, HH and with that of n~ /n + 
ratio somehow leads to a super soft symmetry energy 
x = 1 at suprasaturation densities 0. In the following 
x = is used if not addressed. 

Once the momentum dependence of the single-particle 
potential is obtained, the effective mass of a nucleon with 
isospin t in the nuclear matter medium can be calculated 
from 
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FIG. 1: (color online) Momentum dependence of the single- 
particle potential, the occupation probability and the effective 
mass at p = po and T = 10 MeV ((a), (d), (g)), p = po and 
T = 50 MeV ((b), (e), (h)), and p = 2p andT = 10 MeV ((c), 
(f), (i)) from the MDI interaction. Results from symmetric 
(S — 0) nuclear matter and neutron-rich (S = 0.5) nuclear 
matter are compared. 



In Fig.[T]thc momentum dependence of the single-particle 
potential, the occupation probability and the effective 
mass are displayed at different densities and tempera- 
tures. It is seen that the momentum-dependence of the 
single-particle potential is weaker at higher temperatures, 
which leads to a larger effective mass with increasing tem- 
perature. A higher temperature also leads to a more 
diffusive Fermi surface. At higher densities, the effec- 
tive mass is smaller and the distribution function be- 
comes less diffusive. The effective mass of neutrons is 
larger than that of protons in neutron-rich nuclear mat- 
ter from the MDI interaction. This leads to a symmetry 
potential U syrn = (U n — U p )/8 that is consistent with 
the energy dependence of the Lane potential constrained 
by the nucleon-nucleus scattering experimental data as 
discussed in Ref. [32j]. The relation between the shear 
viscosity and these single-particle properties will be dis- 
cussed in the following. 



III. SHEAR VISCOSITY FROM THE 
RELAXATION TIME APPROACH 

In the present work I study the shear viscosity of 
an isospin asymmetric nuclear matter with uniform and 
static neutron and proton density p n and p p , respectively, 
and temperature T. The static flow field in the nuclear 
system is assumed to be in the z direction and the magni- 
tude changes linearly with x, i.e., u z = cx, u x = u y = 0. 
The equilibrium momentum distribution n is a Fermi- 
Dirac distribution as Eq. in the local frame while it 
is shifted by mu with the flow field in the lab frame. 
Due to collisions, the momentum distribution may differ 
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slightly from the equilibrium one and is written as n and 
6n — n — n° is their difference. 

The shear force between flow layers per unit area is (26j 



F = y^<(p z - mu z )p T -^- > 



^2 d f(p z - mu z )^n T -^^. 



of isospin r moving in the x direction per unit time per 
unit area, where can be written as + (V p [/ T ) x ] 

from Eq. ([5]), and p z —mu z is the momentum transported 
per nucleon. It is easily seen that the contribution of n° 
to the shear force is zero as the integrand is odd in p x . 
Thus, the shear force can be written as 
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In the above p T ^r is the flux, i.e., the number of nucleons 
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To calculate 5n T , we start from the BUU equation using the relaxation time approximation. The isospin-dependent 
BUU equation is written as 
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where the terms 1 — n are from the Pauli blocking effect 
and the degeneracy d — \ takes identical nucleon colli- 
sions into account. In the first-order approximation, the 
contribution of 5n T on the left side can be neglected and 
n T is replaced by n^.. The first term on the left side van- 
ishes as n° is the equilibrium distribution and the system 
is static. By using the relation 
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the second term can be calculated as 
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where U® is the single-particle potential corresponding to 
the equilibrium distribution vP r . The third term on the 
left side can be similarly expressed as 



V r U° ■ V p n° 



Pz dU° 

mc — 

p dp 

So, the left side can be expressed as 
driripx) 
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p dp J 
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dt 

du z p z p x dn° T 
dx p dp 

as c = du z /dx. 
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The right side vanishes if n T = when the detailed balance is satisfied. In the relaxation time approximation only 
Sn T (pi) is kept and the right side can be written as [26j 
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where T T (p\) is the relaxation time, i.e., the average collision time for a nucleon with isospin r and momentum pi, 
which can be expressed as 
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(pi) being the average collision time for a nucleon with isospin r and momentum pi when colliding 



with other nucleons of same (different) isospin, and they can thus be calculated from 



'(Pi 



1 (2tt£ 
2 J (2tt) 3 



p\dp 2 d cos 9 i2d cos 8 



da T T 



in 



P2 



m*(pi) m*(p 2 ) 



x [n° T ip 2 ) - ra?(p 2 K(p'i) ~ n° T (p 2 )n M) + n°M)n°M)] . 



diSS 



(Pi) 



= d 



(27T) 



p\dp 2 d cos 9 1 2 d cos 



d(T T — r 



Pi 



P2 



m%(pi) m*_ T (p 2 ) 



'(2tt)3 y ^ " dft 
x [n°_ r (p 2 ) - n^CpaKCpi) - n°_ T (p 2 )n°_ T (p' 2 ) + n M)n°- T (j>' 2 )] , 

I 



(15) 



(16) 



where 9i 2 is the angel between p\ and p 2 , and is the 
angel between the total momentum p tot = p\ + p 2 = 
Pi + P2 ano - the relative momentum of the final state 
Prei = Pi —p> 2 ■ As for elastic collisions \p\ 1 = \P\ ~P2 1 j 
their magnitude can be calculated from 
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2pip 2 COS 012- 

The magnitude of p'j and pl 2 can then be obtained from 

^Ptot +Plel + 2 PtotPrel COS 8, (19) 
2 \/ Ptot + Plel - 2 PtotPrel COS 9. (20) 
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We choose the isotropic nucleon-nucleon cross sections, 
and in free space they are taken as the parameterized 
forms [33| 

Opp(rm) = 13.73 - 15.04/w + 8.76/v 2 + 68.67w 4 j(21) 
cr np = -70.67 - 18.18/w + 25.26/w 2 + 113.85w,(22) 

where the cross sections are in mb and v is the velocity 
of the projectile nucleon with respect to the fixed target 
nucleon. The above cross sections are the same as used in 
the IBUU transport model calculations. Figure [5] shows 
the pp and np cross sections as functions of the center-of- 
mass energy of the two colliding nucleons from the above 
parameterized form. This parametrization describes very 
well the experimental data for the beam energy from 10 
MeV to 1 GeV [33| , corresponding to the center-of-mass 
energy ^fs from 1883 MeV to 2325 MeV. It is seen that 
for the most probable energies the cross section is smaller 
for pp collisions than for np collisions. 

As used in the IBUU transport model, by assuming all 
the matrix elements of the nucleon-nucleon interaction 
are the same in free space and in the medium [34j , the in- 
medium nucleon-nucleon cross sections are modified by 
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FIG. 2: (color online) pp and np cross sections in free space 
as functions of the center-of-mass energy of the two colliding 
nucleons from the parameterized form (Eqs. (|21|) and (|22[) ). 



the in- medium effective masses in the following form 31 1 

(23) 



^medium „ I ^NN 
a NN = a NN 

' 1-tNN 



where [Inn (i^nn) ^ s ^ ne free-space (in- medium) reduced 
mass of colliding nucleons. 

Eqs. and lead to 



5n T {p) = T T {p) 



du z p z p x dn° T 
dx p dp 



(24) 



As F = — ?7^f , the shear viscosity can be calculated 
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from Eqs. © and (J24J) as 
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where p = yPx + P y + P 2 Z an d n * is the local momentum 
distribution. The second equality sign comes from that 
the shear viscosity is independent of the magnitude of 
the flow and u z = is chosen. 



IV. RESULTS AND DISCUSSIONS 
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The momentum dependence of the total relaxation 
time for a nucleon and that for the nucleon to collide 
with other ones of same or different isospin in symmetric 
and neutron-rich nuclear matter are displayed in Fig. [3] 
A smaller relaxation time means the nucleon on aver- 
age experiences more frequent collisions. A constant 
cross section would make the relaxation time decrease 
with increasing momentum, as higher-momentum nucle- 
ons are more likely to collide with others. Using the 
energy-dependent free-space nucleon-nuclcon cross sec- 
tions, there are peaks around p = 500 MeV, corre- 
sponding to the minimum values of the free cross sec- 
tions as shown in Fig. [2] Due to the smaller effective 
masses of nuclcons in the nuclear medium, which leads 
to smaller in-medium nucleon-nucleon cross sections, the 
relaxation times are larger. In symmetric nuclear mat- 
ter, r r , r* ame and t^ 1 ^ are the same for nucleons of 
different isospins. In neutron-rich nuclear matter, r* ame 
are smaller while T^ ame are larger compared with that in 
symmetric nuclear matter, due to larger chances for nn 
collisions while smaller chances for pp collisions. In addi- 
tion, becomes larger while T p i " L ^ becomes smaller, 
due to a smaller number of protons and a larger num- 
ber of neutrons to collide with, respectively. It is seen 
that there are always more chances for collisions between 
nucleons with different isospins than those with same 
isospin, mainly because of the larger degeneracy factor 
for r di f f than r same as in Eqs. C£5]) and (HH) and the 
larger cross section of <t t ,-t over er r , r . The total relax- 
ation time is dominated by the smaller one, i.e., r dl ^ . 
Thus t„ > t p is found in neutron-rich nuclear matter, 
and they are both of the magnitude only a few fm/c at 
p = po and T = 50 MeV. 

In Fig. 2] the effects of density and temperature on the 
relaxation time arc displayed. It is seen that compared 
with the case at T = 50 MeV, the relaxation time is gen- 
erally larger at T = 10 MeV, due to the fact that the nu- 
cleons are less energetic at lower temperatures. Further- 
more, there appears a peak around p = 270 MeV indicat- 
ing a stronger Pauli blocking effect near the Fermi sur- 
face. This can be understood as nucleons near the Fermi 
surface are more likely to collide with those below the 



FIG. 3: (color online) Momentum dependence of the total re- 
laxation time and the relaxation time for the nucleon to col- 
lide with other ones of same or different isospin in symmetric 
(S = 0) and neutron-rich (8 = 0.5) nuclear matter. 



Fermi surface due to larger cross sections at lower center- 
of-mass energy, but these collisions are largely blocked at 
T = 10 MeV. At p = 2po the peaks move to a higher mo- 
mentum and are even stronger, as the Fermi momentum 
is higher and the Pauli blocking effect is stronger due to 
the less diffusive distribution function as shown in Fig. [1] 
Again the relaxation times are larger from smaller in- 
medium cross sections, and at higher densities and lower 
temperatures this effect is larger due to smaller nucleon 
effective masses. 
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FIG. 4: (color online) Momentum dependence of the total 
relaxation time at different densities and temperatures in 
symmetric nuclear matter from free-space cross sections (left 
panel) and in-medium cross sections (right panel). 
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The specific viscosity, i.e, the ratio of the shear viscos- 
ity over the entropy density, is shown in upper panels of 
Fig.[5]as a function of the temperature for different densi- 
ties (panel (a)) and as a function of the density for differ- 
ent temperatures (panel (b)) in symmetric and neutron- 
rich nuclear matter from free-space nucleon-nucleon cross 
sections. The lower bound rj/s ~ h/Air obtained by the 
Anti de Sitter / conformal field theory (AdS/CFT) cor- 
respondence |35( is also shown by dotted lines for ref- 
erence. The specific viscosity decreases with increasing 
temperature, as a result of more frequent collisions and 
weaker Pauli blocking effect at higher temperatures. Al- 
though the effective mass decreases with increasing den- 
sity, which leads to a larger flux between flow layers as in 
Eq. (|7[l. and a larger relative velocity between nucleons in 
Eqs. (fT5")) and p^|) . the Pauli blocking effect is stronger at 
higher densities. The combined effects lead to an increas- 
ing trend of the specific viscosity with increasing density 
at lower temperatures, but a slightly decreasing trend 
with increasing density at higher temperatures when the 
Pauli blocking effect is much weaker. Furthermore, it is 
interesting to seen that the specific viscosity is larger in 
neutron-rich nuclear matter than in symmetric nuclear 
matter at all the densities and temperatures, especially 
at lower temperatures. This is because that the neutron 
relaxation time is larger than proton in neutron-rich nu- 
clear matter and it contributes more to the total shear 
viscosity according to Eq. (|25|) due to a less diffusive neu- 
tron momentum distribution function. At lower temper- 
atures, this effect is larger because of the larger difference 
between c TiT and <7 t ,-t, stronger Pauli blocking effect on 
neutrons, and less diffusive momentum distribution of 
neutrons as in Eq. (|2~51) . 



The density and temperature dependence of the spe- 
cific viscosity from in-medium cross sections are shown 
in the lower panels of Fig. [5l Compared with the 
results from free-space nucleon-nucleon cross sections, 
the specific viscosity is much larger at lower tempera- 
tures and higher densities, when the in-medium cross 
sections arc smaller from smaller in-medium effective 
masses. The isospin effect on the specific viscosity is 
somehow smaller, compared with the results from free- 
space nucleon-nucleon cross sections, due to the isospin- 
dependent modification on the in-medium cross sections. 
As the effective mass of neutrons is larger than that of 
protons in neutron-rich nuclear matter, the isospin ef- 
fect on t„ is smaller from in-medium cross sections than 
that from free-space cross sections, as shown in panel 
(c) of Fig. El At T = 50 McV, the specific viscosity is 
only about 4^5^, which is already close to the value 
of QGP extracted from the transverse momentum spec- 
trum and elliptic flow using the viscous hydrodynamic 
model (36] . The specific viscosity obtained in this way is 
similar to that extracted from the BUU model calcula- 
tion using the Green-Kubo formulas [22[ , which is around 
20 - 30-^ at lower energies and reduces to as low as 6-P- 
at higher energies. The density and temperature depen- 
dence of the specific viscosity at lower temperatures are 
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FIG. 5: (color online) Specific viscosity as a function of the 
temperature at p = po and p = 2po (left panels) and that 
as a function of the density at T = 10 MeV and T — 50 
MeV (right panels) for symmetric (S = 0) and neutron-rich 
(S = 0.5) nuclear matter from free-space cross sections (upper 
panels) and in-medium cross sections (lower panels). The 
lower bound of the specific viscosity is also shown by dotted 
lines for reference. 



also similar to those in Refs. fill Il9j. 
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FIG. 6: (color online) Shear viscosity r\ and specific viscosity 
ri/s as functions of isospin asymmetry S 2 for different densities 
and temperatures. The solid lines are from the linear fit. 



It will be interesting to study in detail how the shear 
viscosity rj and the specific viscosity rj/s change with the 
isospin asymmetry 6, which is shown in Fig. [5] It is 
seen that both r\ and rj/s roughly satisfy the parabolic 
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approximation 

V(P,T,5) « 7 ? (p,T, ( 5 = 0)+7 ?swm (pT) ( 5 2 , (26) 

(2) ( P , t, *) « (2) (p, r, 5 = o) + (2) ( p t)^7) 

\ S / \ S ' \ S / sym 

where the coefficients rj sym (p_T) and (f) SJ/m (p.T 1 ) arc 

comparable to rj(p,T 7 S = 0) and (^) (p,T,6 = 0), re- 
spectively. The large isospin effect is also observed in 
results from the Brueckner theory [23|. The parabolic 
approximation is good even for very large isospin asym- 
metry, and it seems to be better at lower densities or 
higher temperatures. 

V. CONCLUSIONS 

In this paper I discussed the specific viscosity from 
the relaxation time approach by using the MDI inter- 
action and the nucleon-nucleon cross sections from the 
experimental data modified by the in-medium effective 
masses as used in the IBUU transport model calculations. 
In neutron-rich nuclear matter, the relaxation time of 
neutrons increases while that of protons decreases, com- 
pared with that in symmetric nuclear matter. The spe- 
cific viscosity decreases with increasing temperature be- 
cause of more frequent collisions and weaker Pauli block- 
ing effect. At lower temperatures, the specific viscosity 
increases with increasing density due to increasing ef- 
fect of Pauli blocking while at higher temperatures it 
slightly decreases with increasing density due the smaller 
in-medium effective masses. The specific viscosity in- 
creases with increasing isospin asymmetry, mainly from 
the larger relaxation time of neutrons in neutron-rich nu- 
clear matter. Both the shear viscosity and specific vis- 
cosity roughly follow the parabolic approximation with 



respect to the isospin asymmetry. 

As in this frame work, the effect of the interaction on 
the specific viscosity only comes from the effective mass, 
i.e., the momentum dependence of the single-particle po- 
tential, the value of x, which determines the density de- 
pendence of the symmetry energy, will not affect the re- 
sults. In reality, the transitive matrix T T>T i will also be 
modified in the medium, especially at higher densities. In 
a more realistic Brueckner-Hartree-Fock calculation [33] , 
the in-medium cross section is calculated by consistently 
taking account of the two-body and three-body nucleon 
forces. It was found that the in-medium cross section is 
more isotropic compared with the cross section in free 
space as the forward and backward peaks are largely re- 
duced, and including the three-body force would further 
reduce the nucleon-nucleon cross section. An isotropic 
cross section generally leads to a smaller viscosity than 
an anisotropic one with forward and backward peaks 
as the effective transport cross section is larger in the 
former case, while including the contribution from the 
three-body force seems to increase the viscosity. In addi- 
tion, the system is assumed to consist of uniform nuclear 
matter in the present study. At higher densities and/or 
temperatures, inelastic nucleon-nucleon collisions such as 
JVJV — > NA become important, and A resonances and 
pions will be abundantly produced. At lower densities 
and/or temperatures, liquid-gas phase transition will oc- 
cur and clusters will be formed. All these will affect the 
viscosity of the system. 
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